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Abstract 

We show that the Random Energy Model has interesting rejuve- 
nation properties in its frozen phase. Different 'susceptibihties' to 
temperature changes, for the free-energy and for other ('magnetic') 
observables, can be computed exactly. These susceptibilities diverge 
at the transition temperature, as (1 — T/Tc)~^ for the free-energy. 



A small temperature change in the low temperature phase of spin-glasses 
is able to 'rejuvenate' an already aged system |1|, ^ |^. More precisely, the 
a.c. susceptibility of a system aged for a very long time at Ti < Tc (where 
is the spin-glass phase transition temperature) and then suddenly cooled at 
Ti — AT is, provided AT is not too small, very close the the susceptibility 
of a young system directly cooled from high temperatures to Ti — AT. This 
'fragility' to temperature changes has been interpreted early on as a signature 
of 'temperature chaos', that is the fact that the equilibrium states of a disor- 
dered system are very different for different temperatures: beyond a certain 
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length scale £at (which diverges for small AT), the thermodynamical states 
become uncorrelated and the overlap between them tends to zero for large 
system sizes. Such an effect was conjectured in the context of the droplet 
model of spin-glasses, based on scaling arguments |^, and supported by 
Migdal-Kadanoff renormalisation group calculations [||. A similar effect is 
also predicted for pinned interfaces and has recently been checked in 
careful large scale simulations and analytical arguments in the case of the 
1 + 1 directed polymer in random media |[Tl|. However, temperature chaos 
has been recently disproven in the mean- field SK model and not been 
found either in numerical simulations of the 3D Edwards- Anderson spin-glass 
model [P, |l^. This might be due to the fact that the length £at (if it exists) 
involves a large numerical prefactor and is therefore larger than numerically 
accessible sizes. It could also be that although the length scale i^r plays 
a relevant role in the overlap between the states at Ti and Ti — AT, this 
overlap does tend very slowly to zero for large sizes, as in Finally, the 
'temperature chaos' interpretation of the rejuvenation effect has to be com- 
patible with the simultaneous memory that one observes experimentally. A 



scenario for this was recently proposed in |12 



Another line of thought to explain rejuvenation and memory in spin- 
Lower 



glasses is based on 'hierarchical' energy landscape pictures [|13], |4 



temperatures reveal finer details (rejuvenation) while large scale jumps are 
frozen out (memory) [|15|, |^. Numerical simulations of the dynamics on 
Parisi's tree [jl6| or in the Sinai potential [jl^ indeed confirm that these effects 
exist in the absence of 'true' chaos. Even a simple two-level system does 
actually lead to some rejuvenation when the temperature is of the order of the 
energy difference between the levels - simply because the relative Boltzmann 
weight changes when the temperature is changed. In this paper, we want to 
study in details the temperature rejuvenation effect in the glassy phase of the 
Random Energy Model (rem) , for which a number of exact results are known 
[0. This low temperature phase is described, in the replica language, by 
a 'one-step' replica symmetry breaking scheme [|19|. This model is expected 
to be in the same universality class (with possibly minor corrections) as a 
large number of other models pO[, such as the p-spin models (advocated to 
be good models for glasses), the Bernasconi model, the directed polymer (or 
random manifold) in high dimensions 12^, |2^, the unbinding transition of 
heteropolymers |2^ and the problem of a single particle in a logarithmically 
correlated random potential |p3]. Interestingly, a dynamical version of the 



2 



REM naturally leads to aging dynamics |T^, |25|, ^ . 

Here, we want to compute exactly the 'rejuvenation susceptibility' of 
different observables to small temperature changes. Even if the model is not 
'chaotic', it reveals a number of interesting features that may be relevant 
to the present discussion. For example, the 'susceptibility' to temperature 
changes diverges when the temperature tends to the glass temperature of the 
model. Experimental consequences are discussed in the conclusion. 

It is well known that the low temperature phase of the rem is equivalent 



to that of the 'trap' model [^, where M energy states e^, i = 1, ■■■M are 



chosen with an exponential probability distribution: 

P{e) = ^eM--^)- (1) 

Note that e is chosen to be negative. The partition function for this model 
is simply Z(T) = J^iti^i with Zi = exp(|e|i/T). This model undergoes a 
phase transition at T = Tc, where the partition function 'localizes' on a few 
states. More precisely, in the limit M oo, the Boltzmann weights of a 
finite number of states add up to a finite fraction of the partition function 



for T < Tc [27, P8|. Aging is the dynamical counterpart of this localization 



effect: most of the elapsed time is spent by the system in the deepest available 
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As a first definition of the susceptibility to temperature changes, we study, 
following Fisher and Huse the correlation of the free-energy fluctuations 
for two different temperatures. More precisely, we write: 



where ii stands for logZ(Ti) and the overline means that we average over 
the distribution of the energies e. When Ti = T2(l + e) we expect that: 

CF{Ti,T2) = l-KFe\ (£->0), (3) 

where tip defines the susceptibility to temperature changes. The calculation 



of this quantity starts with Derrida's representation of log Z |jT8 
logZ= r-^^ exp(-t)-exp(-tZ) ^ ^.^ T t^-i (exp(-t) - exp(-tZ)) 

JO t h^0+ Jo 

(4) 
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where b has been introduced to ensure convergence of intermediate calcula- 
tion steps. The average over e then involves: 



exp -tZ = (exp -tz)^ = (1 - (1 - exp -tz)) 



M 



(5) 



For large M, only the vicinity oft — will therefore be of importance. Using 
the fact that the random variables z are distributed with a power law-tail 
^^-i-Ai ^j^j^ ^ _ j'/j'^^ one finds that: 



1 - exp -tz ~t_o r(l - ii)t'^. 



(6) 



Using the last result in the previous two equations finally leads to the fol- 
lowing result for the free-energy (in units of Tg) : 



- /X logZ = -logM - 7(1 -ii)- logr(l - /x), 



(7) 



where 7 is Euler's constant. The first term comes from the fact that when M 
is large, the smallest energy drawn from an exponential distribution behaves 
as — Tc log M, plus order one (random) corrections. Therefore, all the fluctu- 
ations involved in the calculation of Cir will be of order 1. The calculation of 
the average of the product of log Z for two different temperatures is a little 
more involved. At an intermediate level of the computation, one finds, to 
order 6°: 



ili2 



r\b) - m 



r(6//ii 



-|- a- 



+ 1 



b^f3J^{l22, 1/a) 



r^^(i-/^i) 



with a — 1 + e, (5 — a)/ and 



a)^ [ — log 
Jo V 



1 + 



r(i-/.) 



f 

Jo 



du 



r'"^(i-/^2) 



exp(— m) — exp(— m"^) 



u 



(8) 



(9) 



Expanding the previous result to order and rearranging the terms finally 
leads to kf{iii), which is too long to write explicitly here. Its dependence in 
III is given in Fig. 1. Of main interest is its behaviour for small temperatures 
and close to the transition point. For /ii ^ 0, one finds that kf goes to zero 
as: 



hf{i^i) ~ 0.905558.. ./X?. 



(10) 
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Figure 1: Plot of the rejuvenation susceptibilities kq (plain line) and Kp (dotted 
line) in semi-log scale, as a function of = T/Tc. Note that both quantities vanish 
at T = and diverge at T = Tc, with different exponents. 



For small temperatures, only the ground state contributes to the free-energy 
both for Ti and T2. Hence, one indeed expects the sample to sample fluctu- 
ations to be strongly correlated. For /^i — > 1, on the contrary, one finds that 
kf diverges: 

3(41og2-l) 0.538802... 



«;f(a*i) 



(11) 



[I -Hi) (l-/ii) 
indicating that close to the 'delocalisation' transition, the system tends to 
occupy rather different states when the temperature is slightly changed. Fi- 
nally, the correlation of free-energy fluctuations between Ti > Tc and T2 < T^. 
vanish as the size of the system tend to inflnity. 

Therefore, except right at the transition, there is no 'strong chaos' in the 
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REM, which would imply that the susceptibility Kp diverges with the size 
of the system. However, small temperature changes do lead to noticeable 
changes in the physical observables. A perhaps more direct way to see this is 
to assign to each state i a certain observable Oj, independent the energy of 
this state. This can be for example, the position of the state in space if the 
model describes (for example) the metastable states of a pinned interface, or 
the magnetisation of a state for spin system. We will assume for simplicity 
that Ylii Oi = 0. For a given set of random energies e^, the thermodynamical 
value of the observable is: 

1 

{0) = -J:^^0^ (12) 

1=1 

The average over disorder of this quantity is always zero. However, a funda- 
mental difference arises between the case T > and T < T^. In the former 
case and for large systems, (O) tends to zero. This is related to the fact that 
the partition function is more or less evenly spread out on all states. On the 
contrary, for T < Tc, (O) for a given system is a finite random quantity. Its 
variance is (for large M) given by: 

i 

Take for example the case where the observable Oi = Xi is the position of 
a particle in a box of size 2L, Xi = L{1 — 2i/M). In this case, the typical 
average position of the particle is (x) ~ y/1 — fiL which scales with the total 
size of the system, compared to (x) ~ L'' with C < 1 for T > Tc. Quenching 
the temperature from above therefore induces a complete rearrangement 
of the equilibrium properties of the system which occurs in a slow, aging way. 

Now, let us how this observable changes when the temperature is slightly 
changed within the glass phase. The calculation proceeds much as above, or 



as in For small temperature shifts e — 0, we find: 



{{Oh-{Ohf = Ko{f.,)e'[j:0^j, (14) 

where the rejuvenation susceptibility is given by a very lengthy expression. 
This quantity is plotted in Fig. 1. Its small temperature behaviour is given 
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by: 

f^oif^i) ~ ^^T^ /^i = 1.21497802.. ./ii, (15) 

io 

where its divergence for /xi — > 1 is given by: 

-Oil^l) - ^^^y (16) 

Therefore, we again find a divergence of the rejuvenation susceptibihty near 
Tc. For T > Tc, we find that kq vanishes as a (T dependent) power of M. 

We then find a very interesting situation: Eq. (|T^ tells us that when 
the temperature is slighlty changed, the mean position of the particle (say) 
has to evolve by an amount of order ^/KeL that is proportional to the size 
of the system (and with a diverging amplitude when T Tc). In this sense, 
rejuvenation is strong since a small temperature change will induce a rather 
large response of the system. However, since the number of states occupied by 
the particle remains finite in the whole glassy phase T < Tc, the probability 
to find the system in the same state at Ti and T2 remains finite in the limit 
of large systems This probability is directly related to kq, and reads: 



"12 = 1 ^ i^o£ ■ (17) 

(Note that Pu = 1 — /ii as it should [^, ^). In the p-spin glass model 
where two states are generically orthogonal, the resulting two temperature 
overlap function is therefore given by: 



P{q, Ti, T2) = (1 - Pi2)S{q) + P12 6{q - q{T,,T2)). (18) 

In summary, we have shown that the Random Energy Model has inter- 
esting rejuvenation properties in its frozen phase. Different 'susceptibilities' 
to temperature changes, for the free-energy and for other (magnetic) ob- 
servables, can be computed exactly. These susceptibilities diverge at the 
transition temperature with different exponents. Since the rem seems to be 
relevant to many physical situations, the mechanism found here is probably 
of broad interest. However, the coexistence of rejuvenation and memory seen 
in the spin-glass experiments cannot be accounted for by the simplest version 
of the REM, because the evolution at T2 will have a significant influence on 
the properties measured at Ti after reheating. One can generalize the rem 
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along the lines of I^O] by allowing a hierarchy of phase transitions Tc,„. As 
argued in 1 1^ and numerically demonstrated in |]16| , each crossing of a phase 



transition T^^n induces a strong rejuvenation signal (much as calculated here) 
with a slow aging dynamics for T < T^^n and a fast return to equilibrium for 
T > Tc^n, accounting for the memory effect. 
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